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Abstract 



A non-linear device R is connected in series with an admittance Y(uS) and a random noise voltage v(t) 
(e.g., arising from thermal agitation) impressed across the combination. Assume that one knows (a) the 
current-voltage function of R, and (b) the admittance Y(u). Then we can obtain approximate statistical 
information about the voltage Vi(t) across the non-linear device: explicit formulas depending only on (a) 
and (b) can be given for the moments of all orders of «i(i), and similarly for its frequency spectrum. 

1. A non-linear device R is connected in series with an admittance Y(u>) and a random noise voltage v(t) 
(e.g., thermal agitation) impressed across the combination. 



Assume that one knows: (a) the current-voltage function of R and (b) the admittance of Y(lu), Then one 
can obtain approximate statistical information about the voltage v\(t) across the non-linear device: explicit 
formulas depending only on (a) and (b) can be given for the moments of all orders of vi(t), and similarly for its 
frequency spectrum. 

2. There are two ideas involved. The first is to express vi(t) in terms of v(i), assuming (a) and (b) known. 
The second idea is then to make use of the random nature of v(t) to get the statistical information about v±(t). 

The first idea employs an operator-series expansion for vi(t), and the second employs known averaging 
processes on Brownian functions. 

3. In order to make matters definite we shall assume throughout the sequel that if the voltage across R is 
Vi(t) the current through it is V\(t) + •(vi(t)) 2 . 

4. Let us denote by A'(t) the indicial admittance corresponding to the frequency admittance Y(uj). Then 
it is well known that the current through Y(to) is 
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Figure 1. 




Since current through R = current through Y(u), 




Collecting terms in V\ (t) : 
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5. How we some to the first basic step of the paper, that of solving for v\(t) in terms of v(t). To do this we 
assume 1 that 

/oo />oo />oo 

Qi(t - r)v(r)dr + / / Q 2 (t - n,t - r 2 )v(ri)i; (r 2 )dr 1 dr2 
OO J —OO J —OO 

I' oo /*oo /*oo 

+ / / / Q3(*-Ti,*-T2,t-T 3 )t;(ri)u(T2)v(T3)dridr2dT3 (2) 



oo ./ — oo ./ — oo 



it is clear that we shall have secured the desired solution for V\(t) as soon as we have found Q\, Q 2 , Q3, ■ ■ ■ 

6. Now substitute expression (2) for v\(t) in Eq.(l). We shall then equate linear part with linear part, 
quadratic part with quadratic part, etc. 

7. First for the first degree terms in Eq.(l), we have 



Qi(t- T)v(r)dT + / A'(t-a)da / Qi{a - t)v{t)(1t = I A'(t - t)v(t)(It (3) 

1 J— OO J— OO J— OO 

It is clear that (3) will be satisfied if 

/oo 
A'(t-a)daQ 1 (a-r)=A'(t-T) (4) 
-OO 

Suppose Qi is the Fourier Transform of q\\ then recalling that Y(cu) is the Fourier transform if A f (t), i.e., 

/oo 
(?1 (w)e^ t ^ (5) 
-OO 

/OO 
Yiujy^du (6) 
-oo 

Eq. (4) leads to 
Hence, 



gi(w)(l + 27ry(w)) = Y(u) 



thereby determining Qi(t): 

8. Equating the second degree terms in Eq.(l) to each other we have 
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/oo />oo / />oo 

/ Q2{t-Ti,t-T 2 )v(T 1 )v(T2)dT 1 dT 2 + »( Qi(t - T)v{r)dn 
-co J —oo \J — oo 

/OO nOO nOO 

A'(t-r)dr / Q 2 (t - n,r - T 2 )v(T 1 )v(T 2 )dT 1 dT 2 = 0 (8) 
-oo J —oo J —oo 

Eq. (8) will hold if 

/OO 
A'(t-T)Q 2 (T + T 1 ,T + T 2 )dT (9) 
-oo 

Suppose 

/OO pOO 
/ e^ 1+ ^)g 2 (wi,W2)^ lC L;2 (10) 
-oo J — oo 

i.e. Q 2 (ti,t 2 ) is the double Fourier transform of 52(^1,^2). Recalling Eqs. (5) and (6), Eq. (10) leads to 
q 2 (u>i, ui 2 ) + •qi{coi)qi{co 2 ) + 2itY(lji + u) 2 ) q 2 {u)\, ui 2 ) = 0 



1 In general one will also include a constant term Qo 
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Hence 

= »giK)g 2 (^ 2 ) «y(u)i)y(o)2) 

921 2J 2^>l + w 2 ) + l (1 + 2kY(w 1 ))(1 + 2ttY(uj 2 )){1 + 2t<Y{lo 1 +u J2 )) 1 J 

Note that (72 contains the factor •, indicating the effect of the non-linearity of the device R. 
9. In the same way one computes qa(ooi,uj 2 , W3): 



27rgi(wi)g 2 (a;2 + o; 3 

93(^1,^2,^3) - 



1 + 2itY(uj 1 +oj 2 + w3) 

2i 2 y(^i) Y(w 2 )i> 3 ) 



(12) 



1 + 2^^) (l + 27ry(w 2 ))(l + 27rr(cj 3 )(l + 27rr(a; 2 +a; 3 ) 1 + 27ry(wi + W2 + w 3 ) 
where the <7„'s are defined in analogy with Eqs. (5) and (10), as multiple Fourier transforms of the Q n i.e.: 

Q„(ti,...,t„) = J ... J c 3 ^ ulkTk q n {uj u ...,Lj n )duj 1 ...duj n (13) 

Note that (73 contains » 2 . Similarly g 4 will contain » 3 , and so on. Thus taking higher powers of • into account 
is equivalent to going out farther in the series of q's; this is the characteristic feature of perturbation methods. 

10. At this point we interpose a formula which will be needed soon, it expresses the average of Q n , n even, 
in terms of an average of q n : 

/oo />oo 
••• / Q n (Ti,ri,...,Tn,Tn)dTi...dra 
"OO J —OO 

pOO pOO 

= (27r)5/ ... / 9„(wi,-wi,...,w»,w»)(fwi...(fw» (14) 



11. Eqs. (7), (10), and (12) tell us the first three terms of (2). In other words, we have accomplished (to 
three terms) our first basic task, that of expressing vi(t), the voltage across the non- linear device R in terms of 
the admittance Y(uj) and the voltage v(t) across the entire circuit. 

12. So far we have said nothing about v(t), but we are now ready to make use of the fact that v(t) is a 
random voltage. This will constitute the second step of the paper, and will be accomplished by taking averages 
of the random voltages in accordance with known formulas. In these formulas the average is taken with respect 
to the parameter a which in going from 0 to 1 runs through all Brownian motions. x(t, a) is a properly named 
Brownian motion, whose differential is a random voltage and K(t\, . . . , t n ) is a symmetric function of n variables: 



/>1 />oo />oo 

/ da j ... I K(ti, . . . ,t n )dx(ti,a) . . .dx(t n ,a) = 0 
^0 J —00 J —00 



(15) 



if n is odd, 



/ da ... K(ti, . . . , t n )dx(ti, a) . . . dx(t n , a) 

J 0 J —OG J —CO 

n times 

/OO pOO 
... K(t 1 ,t 1 ,t 2 ,t 2 ,...,tn,tn)dt 1 ...dtn (16) 

-oo J — OO 

s V ' 

|times 

if n is even. 

13. Referring to Eq. (2), we inquire as to the average of vi(t). We apply Eqs. (15) and (16) to the Q's, 
then express the result in terms of the q's by (14), and finally (Eqs. (7), (11), and (12)) expressing the q's in 
terms of the admittance Y(oj) we find that the firs non- vanishing term of the average of «i(t) is 

. r 

y_ 00 (l + 27ry(w))(l + 27ry(-a;))(l + 27ry(a;)) 1 ' 
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14. Similarly the average of (vi(t)) 2 is 

/co 
Qi (uj)qi(—uj)duj 
-co 

{pOO pOO pOO >> 2 

2 / / #2(^1, ^2)^2 (-^1, -0J2)duj\duj2 + / q 2 (uj,-uj)du)> 
J —00 J — 00 J —00 ) 

/CO /*CO 
/ gi(wi)g 3 (wi,a;2,-W3)^irfw2 
CO J —QO 

{pOO pOO pOO 
6 / / / 93(^1, uj 2 , u 3 )q 3 (-LJi, -u> 2 , -uj 3 )dujiduj 2 duj 3 
J —CO J —CO J —CO 

/co />oo />oo ^ 
/ / q 3 (uji,uj2, -u 2 )q 3 (-wi,uj 3 , -uj 3 )du}idu} 2 du} 3 *> + . . . (18) 
-00 J -co J —co J 

15. In the same way the higher moments of vi(t) may be computed, so also ma the moments of v(t) — vi(i), 
the voltage across Y(ui), be computed. 

16. An average of much importance is that of + a); this average is called the auto-correlation 
coefficient, and itd Fourier transform gives the frequency distribution of the square of the voltage. The auto- 
correlation coefficient is the average of 

Ql(t - T)v(T)dT + / Q 2 (t-T 1 ,t-T 2 )v(T 1 ) V (T 2 )dT 1 dT 2 ) + 

-co J — 00 J —00 J 

{pOO /"CO /'OO ^ 

/ Q 2 (t + a - T)v{r)dT + / Q 2 (t + ^ - n, i + <r - T 2 Mri)u(T 2 )dTidT 2 + . . . [ (19) 
J —CO J —QO J —CO ) 

which is equal to 

pOO pOO poo 

2-7T / qi(uj)qi(—uj)e~ ju,a da + 4n 2 / duji / dc*>2<?2(wi, —^1)^2(^2, —^2) 



/CO /"CO 
dwl / dLJ 2 q 2 (uj 1 ,oj 2 )q 2 (-uj 1 ,-LJ 2 )c-^^ 1+UJ ^ + ... (20) 
-00 J -co 

Since the second term does not contain a, it is a constant, representing a DC component. The third term is 
equal to 

/CO pCO 
g-jtTcj / <? 2 (w, u> — u>i)q 2 (— U)i, ui — u)du>i 
-co J —co 

Thus to the frequency spectrum 

27rgi (u})qi(-u) 

present with no rectification (i.e. • = 0), there has been added 



8tt 



/CO 
q 2 (uj 1 ,uj- uj 1 )q 2 (-uj 1 ,uj 1 - Lo)dui 
-co 



17. Critique. The method above, of first solving for the voltage across part of the circuit in terms of the 
entire voltage, and then getting statistical averages, is clearly quite general. The particular application of this 
method given here has two weaknesses however: 

(i) the current-voltage relation of R is over-simplified. 

(ii) in every practical case a filter of finite band-width precedes the rectifier-admittance combination of Fig. 
1. 

The problem in which (i) and (ii) have been taken into account can, and should be, set up and solved. 
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